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Several workers have established that the Larkin domains for two three dimensional nonmetallic 
elastic solids in contact with each other at a disordered interface are enormously large, implying that 
there should be negligible static friction per unit area in the macroscopic solid limit. The present 
work argues that the fluctuations in the heights of the random asperities at the interface that occur 
in the Greenwood- Williamson model can account for static friction. 



It is well known that one must apply a minimum force 
(known as static friction) in order to get two solids, 
which arc in contact, to slide relative to each other. 
Several workers have provided evidence, however, that 
there might be no static friction for nonmetallic crys- 
talline surfaces, incommensurate with each other. Aubry 
showed this for the weak potential limit of one dimen- 
sional Frenkel-Kontorova model[l], and recently He, et. 
al.[2] and Muser and Robbins[3] have shown for weakly 
interacting two dimensional incommensurate interfaces 
that the force of static friction per unit area falls to zero 
as A -1 / 2 in the thermodynamic limit, where A is the 
interface area. Even solids made of identical materials 
are incommensurate if their crystalline axes are rotated 
with respect to each other. Disorder, however, can pin 
contacting solids, just as it pins sliding charge density 
waves[4,5] and vortices in a superconductor [6]. Recently, 
it has been shown that Larkin domains (i.e., domains over 
which the solids are able to distort to accommodate the 
disorder at the interface) for contacting three dimensional 
elastic solids are enormously large compared to typical 
solid sizes [7-9], implying that the force of static friction 
per unit area due to interface disorder should also fall off 
as A~ x / 2 in the thermodynamic limit. In Refs. 2 and 3 
it was proposed that the presence of a submonolayer film 
of mobile molecules at the interface is a requirement for 
the occurrence of static friction between incommensurate 
surfaces, which appears to imply that clean interfaces are 
frictionless. In the present work, it is argued that the 
Greenwood- Williamson model (GW)[10,11] predicts the 
existence of sparsely spaced higher than average asperi- 
ties at the interface which interact more strongly with the 
second surface than with each other. Since these are in 
the strong pinning limit, they will exhibit static friction. 

Here, scaling methods, like those used by Fisher for 
charge density waves (CDW)[12], are used to study static 
friction for non-smooth interfaces. This can be accom- 
plished by minimizing the potential energy of the solid in 
contact with a rigid disordered substrate at z=0 with re- 
spect to the size of a Larkin domain[4], which is expected 
to give qualitatively correct results for the problem of two 
elastic disordered solids in contact. Given that the en- 
ergy density of the elastic solid is given approximately 



by 

(l/2)£'|Vu| 2 + V(t)6(z), (1) 

where E' is an effective Young's modulus and V(r) is the 
potential per unit area of the disordered substrate and 
u(r) is the local displacement of the solid, the energy of 
a single Larkin domain is given by 

E = (l/2)L'L 2 E'a[\V t u'\ 2 /L 2 + 

\du'/dz'\ 2 /L' 2 } - V aL, (2) 

where a is a local length scale (e.g., a lattice constant), 
L is the width and L' is the height of the domain, 
V' t u'| 2 — \du' /dx'\ 2 + \du' /dy'\ 2 , where we assume that 
the local displacement u varies on length scales L and L' 
in the x and y and the z directions, respectively. That is, 
we assume that u(x, y, z) has the form u'(x',y',z'), where 
the function u' varies by an amount of the order of local 
length scales (e.g., a lattice constant) when x', y' and 
z', defined by (x',y',z')=(x/L,y/L,z/L'), each vary by an 
amount of order unity. Here, Vq is a typical value of the 
potential per unit area. When Eq. (2) is minimized with 
respect to L' one finds that L' w L and the energy per 
unit volume within a Larkin domain at the interface is 
given by 

E/L 3 » [(l/2)£'|V'uf - V a]/L 2 , (3) 

(where we use the average value of |V'u'| 2 here) whose 
absolute minimum occurs for infinite L (more correctly L 
comparable to the interface length) for i5'|V'u'| 2 > Vba, 
implying that the static friction per atom decreases as the 
reciprocal of the square root of the surface area. These 
scaling arguments apply equally well to disorder on the 
sub-asperity level due to atomic level defects and to the 
multi- micron length scale level due to asperities. When 
applied at the sub-asperity level, they show that the con- 
tact force between asperities should be proportional to 
the square root of area of contact. On the multi- asperity 
length scale, the atoms are replaced by asperities and the 
defects by contacting asperities. 

If the sliding solid has dimensions normal to the in- 
terface much smaller than those along it, L' can only be 
as large as the thickness. Minimization of Eq. (3) with 
respect to L with L' fixed at the thickness shows that the 
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Larkin length is comparable to the film thickness. Thus, 
the interface consists of many Larkin domains. Since the 
pinning force scales with the number of Larkin domains 
at the interface, the force of static friction per unit area 
is non-zero value, in the large solid limit. 

Fishcr[12] has shown that above the critical dimen- 
sions of 4, charge density waves are not pinned for typi- 
cal impurity strengths, but fluctuations in the impurity 
concentration and strength lead to pinning. (The critical 
dimension for two solids in contact is 3, as seen above.) 
Consider the effect of fluctuations in the defect concen- 
tration for thick solids, by dividing the solid into boxes 
of length L and examining the percentage of blocks at 
the interface of sufficiently large defect concentration to 
be in the "strong pinning" regime, where the substrate 
force on each block dominates over the inter-block elastic 
forces. Consider the parameter h « V /aa 2 , where a is 
the interatomic force constant, a is a lattice spacing and 
Vb is the strength of the potential due to a defect acting 
on an atom. Let n c = c'L 2 be the number of defects 
within a particular block and d (where d > c, where 
c is the average defect concentration for the interface), 
the defect concentration strong enough for the block to 
be considered a strong block. Then the ratio of the in- 
teraction of a typical strong block with the substrate to 
aa 2 is ^(c'L 2 ) 1 / 2 . The interface area surrounding each 
strong block is the total interface area A divided by the 
number of strong blocks at the interface, PA/L 2 , where 
P is the probability of a particular block being a strong 
one. Then L 2 /P is the interface area surrounding each 
block and the typical length L' over which the elastic 
interaction between two strong blocks acts is its square 
root, L' = L/ P 1 / 2 . Then the total elastic energy associ- 
ated with each strong block is the product of the volume 
per strong block=(L') 3 times the elastic energy density, 
which is proportional to |Vw| 2 [which scales as (L')~ 2 ] 
or L' . The criterion for a block to be a strong one is 
hic'L 2 ) 1 / 2 » V, or h » (c'P)- 1 / 2 . Since c'P < 1, 
this violates our previous assumption that h « 1, im- 
plying that such fluctuations cannot result in strong pin- 
ning. There are also fluctuations in the locations of the 
points of contact within the defect potential wells within 
each Larkin domain, it too does not lead to static fric- 
tion^]. 

Defect strength fluctuations, however, can lead to 
static friction, as we shall see. On the multi-asperity 
scale, the surface asperity height distribution is likely to 
produce such fluctuations for sufficiently large surfaces 
even those that are quite smooth. This is the situation 
for the GW model[10,ll,13,14], in which there are elastic 
spherical asperities on a surface with an exponential or 
Gaussian height distribution in contact with a rigid sub- 
strate. The GW model is generally accepted to be a cor- 
rect way to account for Anonton's laws in most cases[ll], 
especially for relatively light loads. Volmer and Natter- 
mann's discussion of static friction[14] is not qualitatively 
different from that of Ref. 10. In the GW model, the to- 



tal contact area is of the order of 
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A c = 2irabN ds<j>(s){s - h), (4) 
Jh 

where <p(s) is the distribution of asperity heights z, where 
s = z/a, where a is a length scale associated with the 
height distribution, and h is the the ratio of the distance 
of the lower part of the bulk part of the sliding solid, 
from the surface in which it is in contact, to a, b is the 
radius of curvature of an asperity and N is the number 
of asperities above a certain size, independent of whether 
they are in contact[10,ll,13,14]. The interaction of a sin- 
gle asperity with the substrate is equal to the product of 
the contact area and a shear strength for the interface. 
Actually, to be consistent with our scaling arguments and 
Refs. 2 and 3 we should assume that the friction force 
on a single asperity is proportional to the square root of 
the asperity contact area, but when this was done, the 
results were not changed qualitatively. 

The energy of the interface consists of two parts. One 
part is the single asperity energy, which consists of the 
interaction energy of an asperity with the substrate plus 
the elastic energy cost necessary for each asperity to 
seek its minimum energy, neglecting its elastic interaction 
with other asperities, which is independent of the asper- 
ity density. The second part is the elastic interaction 
between asperities within the same solid, which depends 
on the asperity density. In order to determine these ener- 
gies, let us model the interaction of the I th asperity with 
the substrate by a spherically symmetric harmonic po- 
tentials of force constant at. Assume that in the absence 
of distortion of the solid, the I th asperity lies a distance 
A? from the center of its potential well. Let be the 
displacement of the £ asperity from its initial position. 
We use the usual elastic Green's function tensor of the 
medium at a distance r from the point at which a force 
is applied at the interface, but for simplicity, we approxi- 
mate it by the simplified form G(r) = (E'r) -1 , where E' 
is Young's modulus [15]. Then the equilibrium conditions 
on the u's are 

u e = {E'a^aeiAe - u e ) + ^\ /•;'/.'.., i '<i ,; A, - u,), 

3 

(5) 

where a is a parameter of the order of the size of the as- 
perity. To lowest order in the interasperity interaction, 
approximate the solution for is 

u, = u? + [1 + (E'ay'ae}- 1 ^(E'Rtj) : o,.A, - u°), 

3 

(6) 

where u° is the zeroth order approximation (i.e., the 
solution to Eq. (5) neglecting the second term on the 
right hand side of the equation). Since the contacting 
asperities are randomly distributed over the interface, we 
can estimate the second term (i.e., the summation over 
j) on the right hand side of Eq. (6) by its root mean 
square (RMS) average which is estimated by integrating 
the square of the summand over the position of the j th 
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asperity, which is in contact with the substrate, over its 
position and multiplying by the density of asperities in 
contact with the substrate p. Since the angular integrals 
only give a factor of order unity, we need only evaluate 
the integral over the magnitude of Rtj, giving an RMS 
value of the sum over R^ 1 of order [pln{Lj a)] 1 / 2 where 
here L is the length of the interface and a is the asperity 
size. For L 1cm and a « 10~ 6 cm, [In^/a)] 1 ^ 2 is of 
order unity. The energy of the system can be written as 

(l/2)^a J -|A J --« J -| 2 + 

3 



(1/2)1?' £ J rf 3 r[|VG(r)|(« j |A j - Uj )|] 2 (7) 

3 J 

It follows from Eqs. (5-7) that the two lowest order non- 
vanishing terms in an expansion of the energy of the sys- 
tem in powers of p 1 ^ 2 are the zeroth and first order ones. 
Since the shearing of the junction at the area of con- 
tact of two asperities involves the motion of two atomic 
planes realtive to each other, the distance over which the 
contact potential varies must be of the order of atomic 
distances. Then, if we denote the width of the asperity 
contact potential well by b, of the order of atomic spac- 
ings, we must choose a typical value for a such that ab is 
of the order of the shear rupture strength of the asperity 
contact junction. Thus, a » E'a. Zeroth order in the 
asperity density in Eq. (7) is of the order of a A 2 , where 
a is a typical value of aj, and A is a typical value of 
Aj. The term linear in p 1 ! 2 is easily shown to be of the 
order of E'a 2 A 2 p 1 / 2 to zeroth order in E'a/a. Since it 
depends on p it represents an interaction energy between 
the asperities. Then, the mean inter-asperity interaction 
is proportional to the square root of the number of con- 
tacting asperities per unit surface area, given by 

poo 

p(h) = (N/A) / ds<P(s). (8), 
Jh 

where A is the total surface area and N is the total num- 
ber of asperities whether in contact with the substrate 
or not. The integral in Eq. (4) divided by the integral 
in Eq. (8), which is proportional to the contact area per 
asperity and the square root of the integral in Eq. (8) 
are plotted as a function of the load, which is given in 
this model by 

/•OO 

F L = (4/3)E'N{b/2) 1 / 2 <7 3 / 2 / dscf>(s)(s - /i) 3/2 , (9) 

Jh 

in Fig. 1. 
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FIG. 1. The curve which is lower at the right is a plot of 
the integral in Eq. (4) divided by the integral in Eq. (8) and 
the curve which is higher on the right is a plot of the square 
root of the integral in Eq. (8) versus the integral in Eq. (9). 
All quantities are dimensionless. 

A Gaussian distribution is assumed here for <fr(s) (i.e., 
<f)(s) = (27r) _1 / 2 e _s / 2 ). Since the square root of Eq. 

(8) drops to zero in the limit of vanishing load, whereas 
Eq. (4) divided by Eq. (8) approaches a nonzero value, 
this implies that the interface will approach the strong 
pinning regime (i.e., the regime in which the asperity- 
substrate interaction dominates over the inter-asperity 
interaction) in the limit of vanishing load. 

Let us now give sample numerical values for some of 
the quantities which occur in the application of the GW 
model to this problem. For example, for a typical plas- 
tic, Young's modulus is 4.0 x 10 9 N/m 2 and the shear rup- 
ture strength is 3.5 x 10 7 iV/m 2 [17]. Following Ref. 11, wc 
choose a — 2.4 x 10 _4 mm and b — 6.6 x 10~ 2 mm, and as- 
sume that there is a density of 4.0 x 10 3 asperities/mm 2 . 
Then by performing the integrals in Eqs. (4), (8) and 

(9) , wc find that for F L /A = 3.98 x lO^N/mm 2 , where 
A is the apparent area of the interface, the total con- 
tact area divided by A is 3.03 x 10~ 5 , and the contact 
area per asperity from the ratio of Eqs. (4) and (8) is 
2.44 x 10~ 5 mm 2 . Also, p(h) 1 / 2 , which is equal to the 
square root of Eq. (8) is 1.11mm" 1 . The mean in- 
terasperity interaction force is approximately equal to 
the derivative of the first order term in p 1 / 2 given above 
Eq. (8) with respect to A or E' 'a 2 pih) 1 / 2 A, where a is 
taken as the square root of the mean contact area per 
asperity divided by n. The mean strength of the force 
acting on an asperity, due to the solid with which it is 
in contact, will be estimated by the product of its con- 
tact area and the shear rupture strength E r . Then, the 
condition for the latter quantity to dominate over the 
asperity-asperity interaction, E r Tra 2 > E' kira 2 p 1 / 2 A or 
E r /E' > 4p(/i) 1 ' 2 A, is easily satisfied by the above cal- 
culated quantities since the right hand side is 4 x 10~ 7 
and the left hand side is 8.75 x 10~ 3 . 

Although for higher loads the system appears to move 
towards the "weak pinning" limit, the latter conclusion is 
most likely incorrect because it does not take into account 
the fact that the distribution of asperity heights contains 
asperities which are much higher than average. These 
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asperities will be compressed much more than a typical 
asperity, making the friction force on them considerably 
larger than average. Since the probability of such asper- 
ities occurring is relatively small, however, they will be 
typically far apart, putting them in the strong pinning 
limit. For example, the probability of the ratio of an 
asperity height to a being greater than a value Hl is 

/>OC 

P(h L )= / ds<j>(s), (10) 

Jh+h L 

whose mean height and hence contact area is propor- 
tional to 

/>oo 

PihL)- 1 ds<j>{s)(s-h). (11) 

Jh+h L 

These two quantities are plotted in Fig. 2. It is seen 
that even for only equal to 1/2 (corresponding to an 
asperity height equal to (l/2)cr), Eq. (11) remains larger 
than the square root of Eq. (10). 

Although it has been argued here that the GW model 
predicts the occurrence of a sufficiently dilute concentra- 
tion of asperities with stronger than average forces acting 
on them due to the second solid to consider the asperities 
to be essentially uncorrelated, this still does not neces- 
sarily guarantee that there will be static friction, since 
it has been argued that even for uncorrelated asperities 
static friction will only occur if the asperities exhibit mul- 
tistability[9,16]. The condition for multistability to occur 
at an interface [9], namely that the force constant due to 
the asperity contact potential be larger than that due to 
the elasticity of the asperity (~ E'a), however, will be 
satisfied, as discussed earlier. 
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FIG. 2. Eq. (11) (the higher curve) and the square root of 
Eq. (10) (the lower curve) are plotted versus the load Eq. (9) 
divided by (4/3)E(b/2) 1/2 a 3/2 . All quantities are dimension- 
less. 

In conclusion, although at first sight it appeared that 
arguments based on Larkin domains indicate that the 
disorder at an interface between two nonmetallic elastic 
solids in contact would not result in static friction, when 
one takes into account the fluctuations in asperity height 



that occur in models such as the GW model, there will 
always exist asperities with greater than average height 
to put them in the "strong pinning regime," resulting in 
static friction. For light loads, even typical height asper- 
ities (as opposed to height fluctuations) can easily be in 
the "strong pinning limit." 
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